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Three heat exchangers exist in modern Stirling engines: a heater, a cooler, and a regenerator. Here a study 
that deals principally with tubular heaters and coolers is carried out. The calculation procedure for the 
oscillating flow heat transfer is presented. Literature sources are studied to find the most suitable 
correlations by comparing them to each other and to the classical turbulent flow correlations encoun¬ 
tered in the literature. The enhancement of heat transfer by means of a few circumferential slots inside 
the tubes and the pressure losses of oscillatory flow are discussed. Non-circular cross-section conduits 
with rectangular and triangular cross-sections are investigated and compared to the smooth circular 
tubes. The increment of the performance of an idealised Stirling engine with slotted heat exchanger 
tubes is compared to the case with smooth ones. The ratio of the gain in the shaft power and pumping 
losses is 2.22. The Carnot efficiency increment is 2.7%. 

© 2012 Elsevier Ltd. All rights reserved. 


1. Introduction 

Convective heat transfer equations for oscillating flow encoun¬ 
tered in Stirling engines have not been published until recent 
decades. The common practice has been that the turbulent corre¬ 
lations, e.g. the Dittus—Boelter equation [1,2], are used since ‘no 
correlations existed up to now to calculate the heat transfer coef¬ 
ficient and friction factor in oscillating flow’ [3], Turbulence takes 
time to form and settle and the picture to date is quite complex 
when heating and pressure variation are added. In the fully 
turbulent and fully laminar range the steady unidirectional 
approximations are fairly good [4j. The laminar range is covered 
well in the latest papers in terms of oscillating convective heat 
transfer. In the literature survey several papers were found about 
heat transfer for the oscillating pipe flow encountered in Stirling 
engine heat exchangers. In this study they are compared to three 
classical unidirectional flow equations (Table 1 ). 

As early as 1929 Richardson and Tyler [5] discovered the so- 
called ‘annular effect’ in an oscillatory flow in a pipe, i.e. the 
maximum velocity in an oscillatory flow occurs near the wall rather 
than at the centre of the pipe, as in the case with unidirectional 
steady flow. This applies especially at higher cycle frequencies [6]. 
Tang and Cheng (1993) [7] (Eq. (1), Table 1) employed multivariate 
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statistical analysis to obtain an experimental correlation equation 
for the cycle-averaged Nusselt number in terms of three similarity 
parameters: Re, Re Ui and the dimensionless fluid displacement A w . 
Their experimental runs were performed for air in a heated pipe for 
the following ranges of parameters: 7 < Re<7000; 7<Re w <180. and 
0.06<A W <2.21. On the basis of an experimental and numerical 
study, Zhao and Cheng (1996) [8] derived Eq. (2), which consists of 
only two similarity parameters, Re u and A 0 . They plotted Nu/A° 85 
as a function of Re w (0<Re w <500), while the dimensionless oscil¬ 
lating amplitude of fluid A 0 variable got the values 8.5,15.3, 20.4, 
and 34.9. De Monte et al. (1996) [9] (3) made a comprehensive 
analytical study of oscillating heat transfer in Stirling machines 
regarding the inside and outside heat transfer of the heater, cooler, 
and regenerator and the effect of heat-exchange effectiveness on 
the performance of a Stirling engine. They applied the empirical 
correlation proposed by Tang and Cheng (1993) [7] and suggested 
that the correlation was valid under the following range of condi¬ 
tions: 11 < Re max <10,995; 1.75<Re u <45. Walther et al. (1998) [6] 
performed a theoretical analysis to evaluate the influence of 
developing flow on the heat transfer associated with laminar 
oscillating pipe flow. Their heat transfer correlation (4) (the integral 
Nusselt number describing the averaged wall heat transfer as 
a function of the dimensionless oscillating frequency, the pipe 
length, and the bulk velocity) was plotted for the range: 1000 < 
Re max < kCmax’’ I 50 < Re,,, < 900; 0.2 < A < 900, where 
Re l S s = 400Vt&Z. 

In this study Eqs. (1)—(4) were first tested by numerical calcu¬ 
lations in a tubular heat exchanger and then compared to three 
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Nomenclature 

Vse 

swept volume of piston 



V 

mean velocity of flow 

A 

cross-sectional area of moving element or flow 

X 

fluid or piston displacement, local coordinate 

A 0 , A r , A 

^dimensionless amplitude of fluid displacement or 

-^max 

amplitude of fluid displacement 


stroke 



Oh o 2 

inside and outside width of rectangular duct 

Creek symbols 

b h b 2 

inside and outside height of rectangular duct 

a 

heat transfer coefficient 

c 

outer leg of triangular duct 

A 

difference 

Q 

Darcy friction factor 

<5 W 

wall thickness 

d 

inner diameter of circular tube 

Cl 

entrance effect correction factor 

dh 

hydraulic diameter of duct or tube 

Vc 

Carnot efficiency 

i 

function in the correlation formula given by Petukhov 

A 

Fanning friction factor 


and Popov 

A w 

conductance of wall 

c 

total conductance 

A 

relative amplitude of fluid displacement, critical length 

L 

tube length 


ratio 

Nu 

Nusselt number 

P 

dynamic viscosity 

P 

shaft power 

P 

density 

Pf 

pumping losses 

Tw 

shear force at duct wall 

Pr 

Prandtl number 

<2> 

heat flow 

Re 

Reynolds number 

<P 

phase angle, drag coefficient induced by sudden 

Qm 

mass flow rate 


changes of flow by geometry 

<Jv 

volume flow rate 

(j) 

oscillation frequency 

RSm 

Re average over half a cycle 



Rcu 

kinetic Reynolds number (= pojd 2 /n ) 

subcripts 


T 

temperature 

i 

inner 

T h 

heat source temperature 

0 

outer 

T l 

coolant temperature 

c 

cooler 

Um 

cross-sectional mean velocity 

E 

heater 

f^max 

amplitude of the cross-sectional mean velocity 

DE 

dead volume 

U 

mean velocity of unidirectional flow 

F 

friction 

V E 

momental expansion volume 

SE 

swept volume 

+E 

dead volume of expansion space 

X 

local 


classical correlations (5)—(7) by Dittus and Boelter [10], Petukhov 
and Popov [11], and Gnielinski [12], respectively, meant for unidi¬ 
rectional turbulent developed flow but applied to the oscillating 
pipe flow conditions. Eq. (7) of Gnielinski [12] covers a large Re 
range: 3000 < Re < 10 6 . 

The reciprocating flow in Stirling engines changes its direction 
constantly. During one cycle the laminar, transitional, and turbulent 
forms of flow modes occur. In this paper a high-frequency (8-Hz) 
low heat source temperature engine is considered. The review of 
these equations is performed in order to find tools to enhance the 


heat transfer of oscillating flow by narrow boundary layers on the 
walls of heater and cooler tubes. In addition to the circular tubes, 
circumferentially slotted and non-circular rectangular and trian¬ 
gular tubes are studied to improve the performance of the Stirling 
engine that is being studied. 

2. Assumptions and input values 

In order to calculate the heat transfer in the heat exchangers the 
following assumptions are made according to de Monte et al. 


Table 1 

Oscillating and classical unidirectional heat transfer equations. 


Nu 


-0.494 + 0.0777- 


A M 

1+A; 


■Re 07 


0.00162-Re° 4 Re° 8 A a 


Xmax _ rr u m _ d Re /-j\ 
~TT ~~ TjX ~~ iReO, 


Nu = 0.02-A° 85 Re° 5s A 0 = 

Nu = -0.494 + 0.0777 ■( - Ar a -Re° 7 - 0.00162-Re° 4 - (4 Re„) as A R = 

VI +A r J 4 L Re„ 

Nu = 5.78 + 0.00918-Re 8969 /I -0 - 367 -f 0.1781— A~ 0703 Re° 526 A = 

| R&max | 2 u R^max 

Nu = 0.023-Re° 8 Pr° 4 

Nu = - (f/2)RePr - f = (3 . 6 41ogRe _ 3 , 2 8)~ 2 

(1 +13.6 f) + (11.7 + 1.8-Pr '/3)+j (pr2/3 _ !) 


Nu 


(f/8) (Re - 1000) Pr 
1 +12.7-If/8)(Pr 2 / 3 - 1) 


3000 < Re < 10 6 


/ = (0.79 TnRe - 1.64) -2 
10 4 < Re < 5-10 6 


( 2 ) 

(3) 

(4) 

(5) 

(6) 


(7) 
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(1996) [9]: 1) the working fluid temperatures in the heater and 
cooler are uniform and time-independent; 2) the ideal gas equation 
applies; 3) the cycle gas density in the heat exchangers is constant; 
4) an ‘alpha’-type engine is considered; 5) the machine is working 
in cyclic steady-state operation, characterised by sinusoidal laws of 
motion of the piston. Input values are presented in the context of 
the calculations. 

The alpha-type engine is one of three main types of Stirling 
engine arrangements. The alpha-type engine uses two pistons that 
mutually compress, expand, and move the working gas between 
hot and cold spaces (Fig. 1). See, for instance, Martini 1983 [13] for 
a more detailed description of the arrangements. The tubular heat 
exchangers in the engine that are studied are realised by counter¬ 
flow configuration. 

3. Piston motion and dimensionless numbers 

The time-averaged velocity in a Stirling engine is zero, or at least 
very low, since in a Stirling engine we do not have a pulsating flow, 
but an oscillating flow, meaning that the stationary component of 
the flow rate is zero. 

The Schmidt theory is one of the isothermal calculation 
methods for Stirling engines [14]. For the alpha-type engine the 
momental expansion volume 14 in Eq. (8) is function of a swept 
volume of the piston Vse and the dead volume of expansion space 
Vde [14]. 


where q> is the crank angle. By substituting V = x-A and A E = Ase = 
■Ade = A, we get 

x E -A = Xn y A (l-cos(wt))+x DE -A (9) 

Differentiating Eq. (9) with respect to time gives the cross- 
sectional mean velocity u m 

u m = ^|^.sin( W t) = ^f^sinrp (10) 

The phase angle <p = ut of the cross-sectional mean velocity 
ranges from 0 to 2tc. The amplitude of the cross-sectional mean 
velocity u max (Eq. (11)) is obtained when <p = tc/2. It is related to the 
amplitude of the fluid displacement x max by the piston [8] (stroke, 
bore: 0.15 m, 0.3 m) and is the distance between the top and the 
bottom dead ends. 


Eq. (10) is used as the hydrodynamic boundary condition for the 
calculation of the heat transfer characteristics in the oscillating pipe 
flow. The mean velocity ti m is illustrated in Fig. 2 as a function of the 
crank angle (where u max is 3.75 m/s). 

By means of u m and u max we can calculate Re m , Re max and q m m 
and q m max . These parameters are encountered in the literature 
correlations (1)—(4) of the oscillating pipe flow. 

Fig. 2 presents the cross-sectional mean velocity for piston. For 
the heater tubes that were studied (helium as working medium: 
150 °C, 120 bar) the corresponding cross-sectional mean velocity 
amplitude u max is 14 m/s and Re max = 47,423 at a frequency of 8 Hz 
because of the different face areas between the piston and the 
bundle of 670 tubes. According to de Monte et al. (1996) [9], 

14 m = —Re ma x (12) 

7T 

The Re (=Re m ) linked to the Eqs. (1) And (3)—(7) is based on the 
cross-sectional mean velocity u m . It is a kind of average over half 
a cycle (Eq. (12)). 

4. Calculation procedure 


V E = (1 - COS rp) + V DE 


( 8 ) 


Eq. (2) by Zhao and Cheng (1990) [8] is used in demonstrating 
the calculation sequence. u m is calculated from the Eq. (10), while w 
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Fig. 2. Cross-sectional mean velocity u m as a function of the piston crank angle. 
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Fig. 3. Nu/4„’ 85 as a function of Re„. 
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and the stroke x max are constants of 50.3 1/s (~8 Hz) and 0.15 m 
respectively. The phase (crank) angle of the piston <p and u m vary 
according to Fig. 2. Re m also varies but Re a has a constant value of 
1022. The thermophysical properties of helium at a pressure of 
120 bar [15] are used to calculate the density p, the dynamic 
viscosity p, the Prandtl number Pr, and the conductivity A that 
appear in the dimensionless numbers. By using the mean velocity 
amplitude of 14 m/s encountered in the heated tubes (u max ), we can 
calculate the amplitude of the fluid displacement in the tubes 
(x max ) from Eq. (11). The dimensionless oscillating amplitude of 
fluid A 0 in the pipes is calculated through the pipe diameter d and 
Xmax- Finally Eq. (2) gives the mean value of 52.6 W/(m 2 K) for Nu. 

According to Fig. 2, we see that u m receives its maximum value 
when the phase angle is 90°. By changing the frequency 0 </< 8 Hz 
(0 < w < 50 1/s) and keeping q> as a constant (7t/2), we can plot 
Nu/Ao 85 as a function of ke w [8], as presented in Fig. 3. 

The plotted Nu/A ° 85 as a function of Re u is well compatible with 
Zhao and Cheng 1996 [8] but it is extended here to the range 
0 < Re u < 1022, while it was originally presented for the 
range ~ 0...500. 

In this paper we will plot the graphs for Nu as a function of Ke w . 


5. Heat transfer 

Pressurised hot water (T H ) on the shell side is used as a heat 
source, as shown in Fig. 1. For the heat transfer calculations 
a constant temperature of 150 °C and a pressure of 120 bar inside 
the heater tubes were selected for the working medium, helium. 
The Nu numbers of correlations (1)—(7) on the helium side are 
plotted in Fig. 4 as a function of Re^. 

Fig. 4 indicates that when the frequency of the oscillation is 
small some of both correlation groups (oscillating and unidirec¬ 
tional turbulent flow) have almost the same slopes and the graphs 
are almost convergent. As the Re w increases the slope of the laminar 
oscillating correlations starts to decrease gradually. 

Eq. (2) from Zhao and Cheng (1996) [8] gives the second largest 
Nu values of the oscillating graphs. It is a conservative estimate of 
the heat transfer rate for a reciprocating flow in a pipe of infinite 
length. Here only two dimensionless parameters, A 0 and Re Ul are 
included. The other oscillation flow Eqs. (1), (3), and (4), which 
include the classical Re number as a third non-dimensional term 
(commonly used in the unidirectional pipe flow) give smaller Nu 
values. 


O Zhao & Cheng □ Tang & Cheng mod. A Walther et al. mod. X Walther et al. 

Xde Monte et al. O Dittus & Boelter +Gnielinski - Petukhov & Popov 
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Fig. 4. Nu numbers of correlations studied as a function of Re w . 
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By comparing the non-dimensional oscillating amplitudes A R 
and Ay in the sources de Monte et al. (1996) [9] and Tang and Cheng 
(1993) [7], it is evident that Re = Re max . The modified (mod) 
diagrams are arrived at by this assumption. Under this circum¬ 
stance the highest oscillating Nu values were plotted from Eq. (1). 
Otherwise, we did not get the positive Nu values from that 
equation. 

We did not find the oscillating heat transfer correlations for the 
turbulent flow. However, we can conclude that unidirectional 
turbulent correlations give slightly higher dimensionless heat 
transfer coefficients at the higher frequencies (Re w numbers) than 
the calculated coefficients from the oscillating equations of Zhao 
and Cheng (2) and the modified Tang and Cheng (1). 

Ahn and Ibrahim (1992) [16] conducted two-dimensional 
analysis of oscillating gas flow in laminar, transition and turbu¬ 
lent regimes inside Stirling engine heat exchangers. Numerical data 
were compared to the experimental results of other investigators 
with good accuracy in laminar and turbulent flow regimes. They 
reviewed several criteria for the critical Reynolds number for 
oscillating pipe flow. The onset of the turbulence and frictional 
losses in an oscillatory turbulent pipe flow is discussed in another 
source, Zhao and Cheng (1996) [17]. The oscillating flow is in the 
turbulent region when A 0 y/ReZ > 961 [17]. At the 8-Hz frequency 
that was studied we get 2966 for such a product (of A 0 = 93 and 
Re w = 1017) and the flow is turbulent. 

6. Heat transfer enhancement 

In this work several traditional and some up-to-date heat 
transfer enhancement methods, e.g. regarding the laminar 
convective heat transfer problem as a conduction problem [18], 
were gone through. However, here only the effect of circumferen¬ 
tial slots in a circular pipe and the enhancement of heat transfer by 
non-circular square and triangular ducts are reported. 



According to Mills (1995) [12], turbulent flows with simply 
defined hydrodynamics at a pipe entrance are seldom encountered 
in engineering practice. Most often there is a sharp 90° edge, 
a bend, or an elbow (Fig. 5). 

The corresponding hydrodynamic entrance lengths vary from 
about 10 to 15 diameters, where no large-scale eddies are present, 
to about 30—40 diameters, where there are large-scale eddies [12]. 
The effect of entrance configuration on the average heat transfer for 
the turbulent tube flow is presented in the Fig. 6. 

It is visible in the Fig. 6 that the open-ended 90° edge and 90° 
elbow give the largest eddies and the most long-range average Nu 
enhancement. The observation that ‘the edges improve the inside 
heat transfer’ in the tubular heat exchanger was studied more 
carefully. 

It is suggested that a few sharp edges inside the pipe could 
produce a greater local Nu number and a narrower laminar 
boundary layer along the entire length compared to the smooth 
pipe. The axial slots (grooves) are used e.g. in the annular ducts of 
electric machines on the stator and/or rotor surfaces in order to 
enhance the heat transfer [19]. 

The heater tube length (170 mm) was divided into four similar 
parts by means of three circumferential slots (Fig. 7). 

The depth and width of one slot opening are assumed to be 
2—3 mm and 7.5 mm, respectively. In this case the effective heat 
transfer length is still 170 mm, even if the total pipe length with 
slots is increased to 192.25 mm. 

A factor and numerical fit was made for the open-ended 90° 
edge entrance effect. Eq. (2) from Zhao and Cheng (1996) [8] was 
multiplied by this factor. The effects in the entrance and after the 
sequential slots in the straight plain tube were assumed to be 
identical. In addition this was compared to the factor ei in Eq. (13) of 
Sukomel et al. (1987) [20] (Fig. 8). 
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Fig. 6. The effect of entrance configuration on the average heat transfer in the turbulent tube flow. 
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OZhao& Cheng Dx EpsL, no slots Ax EpsL (Sukomel et al.) 



x/d 


Fig. 8. Local Nu x enhancement divided by cycle-averaged Nusselt number Nu. 


= 0.022Re° 8 Pr 043 ei 
= 1.38(x/dr 012 


Nu x 


In Fig. 8 we can see the mean heat transfer from Zhao and Cheng 
(2) and the entrance effect without and with circumferential slots. 
The effect of the factor e\ [20] with slots compared to the mean heat 
transfer looks promising. It is selected to account for the change in 
the heat transfer coefficient in the initial section of the tube and 
after the slots. 

The performance of the selected non-circular cross-sections of 
rectangular and triangular tubes was investigated and compared to 
that of the circular tube. Eqs. (14)—(16) for the total conductance G 
were derived for the circular, rectangular, and equilateral triangular 
tubes. 


7. Pressure losses 


The internal flow pressure loss of the smooth and slotted tubes 
was studied and the loss of the smooth tubes of a non-circular 
cross-section (rectangular and triangular) was compared to the 
loss in the circular tube. 

Two different friction factor definitions are in common use in 
the literature [21], the Darcy friction factor Cf and the Fanning 
friction factor A: 


Q 


2t w 1 ^ 
pu 4 


The losses are calculated as follows: 


(17) 


l_ 

C 


1 1 1 

$j7rdj Xyy dj + do a 0 7rd 0 

5w 


(14) 


1 is the pipe length, A w and <5 W the conductance and thickness of the 
wall. For the circular pipe dj and d 0 are the inner and the outer pipe 
diameters and cq and a 0 the corresponding heat transfer coeffi¬ 
cients at the surfaces. 


G 2a i (a 1 +a 2 ) + a : + a 2 + hi + b 2 + 2a 0 (b 1 + b 2 ) 

-^A/V-C-~- 

b i - ai 

For the rectangular duct a\ and a 2 (and bi and b 2 ) are the inside 
(and outside) width and height of the duct. 


I _ 1 1 1 

C o ( 2i5 w \ A w / 5 W \ 3 a 0 c 

v tan30°y <5wV tan30”J 


(16) 


A P . f (18) 

The drag coefficient q> includes all sudden changes in the flow as 
a result of the geometry, dh is the hydraulic diameter. 

Eq. (18) applies all kinds of flow patterns, e.g. laminar and 
turbulent. The drag coefficient A depends on the Reynolds number 
for the flow within the duct. The hydraulic diameter dh allows Eq. 
(18) to be applied to any given cross-section [22], The plot of A is 
presented in the source [22] for tubes of various cross-sections as 
a function of the Reynolds number. 

The friction coefficient Cf of unidirectional flow can be calcu¬ 
lated from the correlation (19) for laminar and from the correlation 
(20) for the turbulent flow [23]. 

C f = ^, Re< 2000 (19) 

1 Re 

C f = 0.079 Re-*, Re < 10 5 (20) 


For the triangular equilateral duct c is the outer leg of the 
triangle and 5 W the wall thickness. 

It was found out that the smooth non-circular square and 
triangular tubes have greater free flow and exterior areas than the 
circular tube with an equivalent hydraulic diameter [4[. Hence, for 
the equivalent heat flow magnitude the required tube lengths of 
square and triangular ducts were 21% and 39% shorter than the 
length of the circular tube. However, the dead volume caused by 
the heat exchanger remains constant. 


The friction coefficient of a cyclically turbulent oscillatory flow 
was measured by Zhao and Cheng in 1996 (Eq. (21)) [17]. 

Q =i(S? +040624 ) ,21 > 

which is valid in the range of 81 < Re w < 540 and 54.4 < A 0 < 113.5. 

The laminar and turbulent drag coefficients for unidirectional 
tube flow were plotted when A = 4Cf as a function of Re u . This curve 
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Fig. 9. The laminar and turbulent drag coefficients for unidirectional tube flow compared to the drag coefficient of a cyclically turbulent oscillatory flow. 


is compared to the friction coefficient graph of the cyclically 
turbulent oscillatory flow in Fig. 9. 

The oscillatory Eq. (21) for the turbulent flow gives slightly 
lower values for the friction coefficient than the unidirectional 
friction Eq. (20) (Fig. 9). In the studied heat exchanger that was 
studied A 0 = 93.3 and Re u = 1016.8 (~ 8 Hz) and the flow is 
turbulent, after Zhao and Cheng [17]. 

The oscillating pressure loss for the circumferentially slotted 
circular tube (3 slots, 2 mm depth) is 8.17 times higher than for the 
smooth tube at a frequency of 8 Hz. The friction loss towards one 
tube is small (~ 1 W), however. 

The pressure losses of unidirectional and oscillating flow for the 
smooth circular and non-circular tubes are presented in the Fig. 10. 
It is visible that unidirectional flow in the tube retains larger losses 
than the oscillating case. The length of the circular tube is 0.17 m 
and the inner diameter is 0.006 m. 

The square and triangular tubes correspondingly have only 
marginally lower losses (as a result of the diminished pipe lengths 
that give equivalent heat transfer) compared to the smooth circular 
tube. 


Friction coefficients and Nusselt numbers for fully developed 
laminar flow are presented for a large variety of tube shapes in the 
source Rosenhow et al. (1985) [24]. In all cases the hydrodynamic 
diameter is used as the characteristic length in the Reynolds and 
Nusselt numbers. The friction CfRe and Nu values of the regular 
polygonal tubes decrease in a gently sloping manner as the number 
of sides decreases. 

Bejan and Lorente (2006) suggest that the round tube shape is 
the best regarding the laminar flow resistance. However, nearly 
round shapes such as hexagonal and square tubes perform almost 
as well after the constructal law. Even if the duct cross-section is 
imperfect — that is, with features such as sharp corners, which 
concentrate fluid friction — its performance is already as good as it 
can be [25], 

8. Performance with slotted tubes 

The thermodynamic performance of an idealised Stirling engine 
was studied by comparing the heat transfer between circum¬ 
ferentially slotted and smooth circular tubes in the heat 
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Fig. 10. Pressure losses in one smooth tube: unidirectional and oscillating flow in the circular tube and the oscillating flow in non-circular square and triangular tubes. 
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exchangers. The identical heater and cooler by the geometry for the 
engine that was studied were roughly dimensioned. The effect of 
the slots was considered from two points of view. 

First, the difference in the shaft power was calculated according 
to Eq. (22) when comparing slotted and smooth tubes in the heater. 

A P = t)-A<P E (22) 


heater (and conductance Gq and G E ), we got a T E increase of 2.5 °C, 
a Tc fall of 2.3 °C, and Arj ~ 0.0272. It means a 2.72% increment in 
the efficiency of the Stirling engine by using these circumferentially 
slotted tubes instead of smooth ones. 

9. Conclusions 


where A P is the difference in the shaft power when using slotted 
and smooth tubes, A$ E the corresponding heat flow difference in 
the heater and r) the efficiency (= 0.22) of the selected Stirling 
engine. In this first study the logarithmic temperature difference 
AT e between helium and water in the heater was supposed to 
remain constant. The heater heat flow <2> E was increased as a result 
of the rise in the conductance G E and the inside heat transfer 
coefficient cq (Eqs. (14) and (23)) while using the circumferential 
slots compared to the smooth case. 

<2> e = G E AT E (23) 


The <P E values were calculated for the slotted and smooth heat 
exchanger tubes used in the Eq. (22). 

The corresponding (oscillating flow) pressure losses were 
calculated by Eqs. (17), (18), (21) and Eq. (24) 

AP F = q V 'Ap (24) 

for a heater consisting of 670 tubes. The drag coefficients <p was 
taken from the source [26]. As a result we got the power ratio 
(Eq. (25)) 


AP _ 1336,8 W 
APp - 603 W 


(25) 


Here the factor ei (Eq. (13)) [20] together with Eq. (2) from Zhao 
and Cheng [8] was used to calculate the (arithmetic) mean heat 
transfer coefficients for the slotted pipes (3 slots, depth 2 mm) and 
the difference A$ E and AP. 

From Eq. (25) we can see, on the basis of previous assumptions, 
that the gain in the shaft power of the Stirling engine compared to 
the increase in the pumping losses is 2.22 when these circum¬ 
ferentially slotted tubes are used. The same ratio is 1.94 for the 
tubes consisting of 3 slots with a depth of 3 mm. 

We also considered the efficiency increment. In this case we 
kept the cooler and heater heat flows 4>c and $ E and the shell-side 
(water) temperatures constant (Plotted pipe = ^smooth pipe). The 
subscript C refers to the cooler. As a result of the slots in the tube 
both Gc and G E increase and the logarithmic ATc and AT E are 
reduced compared to the smooth tube (Eq. (23)). Hence the new 
working fluid temperatures T E and Tc can be solved for the helium 
in the slotted heater and cooler and the Carnot efficiency (Eq. 
(26)) [27], 


Vc = 1 - 


Tc 

T e 


(26) 


The increment in the efficiency of the Stirling engine can be 
calculated as 


At? c 


Vc slotted Vc smooth 
Vc smooth 


(27) 


where ?j c smooth is the solved Carnot efficiency with a smooth heater 
and cooler tubes and j j c slotted is the calculated efficiency with 
slotted tubes. 

By using the heat transfer enhancement by the factor ei [20] 
together with Eq. (2) from Zhao and Cheng [8] to calculate the 
working fluid heat transfer coefficients in the slotted cooler and the 


Earlier fully laminar and turbulent steady unidirectional equa¬ 
tions were used as an approximation for convective heat transfer in 
a Stirling engine. Hitherto the literature has dealt with the laminar 
oscillating correlations. In this study the existing laminar equations 
are compared and extended to the turbulent range. It is noticed that 
some of them give almost the same Nu values as turbulent unidi¬ 
rectional equations. 

The salient parameters are the amplitude of the cross-sectional 
mean velocity, the kinetic Reynolds number, and the dimensionless 
oscillating amplitude. The equations based on these two latest 
parameters give Nu values that are almost as high as the traditional 
equations in the turbulent range. 

Heat transfer enhancement is studied with a few circumferen¬ 
tial slots inside the smooth tubes. The effect of edges in the 
entrance and the sequential slots is assumed to be identical. In 
addition, non-circular rectangular and triangular tubes are inves¬ 
tigated. The square and triangular tubes have greater free flow and 
exterior areas than the circular tube with the same hydraulic 
diameter. Hence the duct lengths are also shorter. However, the 
dead volume in the heat exchanger remains constant. 

The friction coefficient of oscillating flow is lower than in the 
unidirectional approximation. The frictional losses in slotted tubes 
are greater but with square and triangular tubes they are margin¬ 
ally lower than with smooth circular ones. 

The increment of the performance in the idealised Stirling 
engine with slotted heat exchanger tubes in oscillating flow is 
compared to smooth tubes. The ratio of the gain in the shaft power 
and pumping losses is 2.22 compared to smooth tubes. The Carnot 
efficiency increment is 2.72%. 

With these simplified examinations it can be concluded that 
methods (few edges or circumferential slots) that form narrower 
boundary layers along the length of the tube but do not signifi¬ 
cantly increase the dead volume or pressure loss are recommended 
in the design of the heater and cooler of the Stirling engine. 
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